HEAVY HEXAQUARKS IN A CHIRAL CONSTITUENT 

QUARK MODEL 



S. Pepin* and Fl. Stancu 
Universite de Liege, Institut de Physique B.5, Sart Tilman, B-4000 Liege 1, Belgium 

(February 1, 2008) 

Abstract 

We discuss the stability of hexaquark systems of type uuddsQ (Q=c or 
b) within a chiral constituent quark model which successfully describes the 
baryon spectra including the charmed ones. We find these systems highly un- 
stable against strong decays and give a comparison with some of the previous 
literature. 

I. INTRODUCTION 

According to QCD rules, systems such as tetraquarks, pentaquarks or hexaquarks can, 
in principle, exist. Their study is important in disentangling between various QCD inspired 
models. Here we are mainly concerned with non-relativistic models, which simulate the low- 
energy limit of QCD. In these models, the central part of the interquark potential usually 
contains a linear term which describes the QCD confinement and a Coulomb term generated 
by the long-range one-gluon exchange (OGE) interaction. The spin part is usually described 
by the chromo-magnetic part of the one-gluon exchange interaction, analogous to the Fermi- 
Breit interaction of QED g,||. 
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An interest in the constituent quark model has recently been revitalized after recog- 
nition of the role of the spontaneous chiral symmetry of the QCD vacuum. This implies 
that the valence quarks acquire a constituent dynamical mass related to the quark conden- 
sate < qq > and that the Goldstone bosons 7T, K, r\ couple directly to constituents quarks 
||. It has been shown that the hyperfme splitting and especially the correct ordering of 
positive and negative parity states of baryons with u,d and s quarks are produced by the 
short-range part of the Goldstone boson exchange (GBE) interaction |3|,|||, instead of the 
OGE interaction. 

In Ref. 0, we studied the stability of the H-particle, a uuddss system, with J p = + 
and 1=0, in the frame of the chiral constituent quark model of Ref. ||. We found that 
the GBE interaction induces a strong repulsion in the flavour singlet uuddss system with 
J p = + and 1=0, i.e. this system lies 847 MeV above the AA threshold. This implies that 
the model of Ref. || predicts that the H-particle should not exist, in contrast to Jaffe's || 



or many other studies based on conventional one-gluon exchange models [K| . In the model 
used by Jaffe, the chromomagnetic interaction, gave more attraction for the flavour singlet 
state than for two well separated lambda baryons. In Jaffe's picture the H-particle should 
be a compact object, in contrast to the molecular-type deuteron. 



In a recent study by Lichtenberg, Roncaglia and Predazzi [[L1J| , the uuddss system is 
discussed in the context of a diquark model and it is found unstable, in contrast to Jaffe's 
result. The above authors also discuss hexaquarks in which one of the s quarks is replaced 
by a c quark, the H c particle, or by a b quark, the particle. The charmed hexaquark is 
found unstable but the bottom hexaquark is found stable by about 10 MeV with respect to 



the A + Ab threshold. Based on the concept of dynamical hadron supersymmetry, Ref. JTT 
predicts the exotic masses from the ones of ordinary mesons and baryons as input, with no 
free parameters. On the other hand, more sophisticated calculations within a constituent 
quark model with chromomagnetic interaction give both H C (1=0, J=3) and ifj,(I=0, J=2 
or 3) stable by 7.7 MeV up to 13.8 MeV ||12| . In fact, from general grounds |I3"| , [H| one 



expects that the stability of multiquark systems should increase with the mass asymetry 



of the constituent quarks and this was tested for tetraquarks systems, in a conventional 
constituent quark model with chromomagnetic interaction [p~5] p~T[] . 

Here we focus our attention on hexaquarks and in particular on the uuddsQ (Q=c or 
b) system, which results as a promising candidate from the above models. We study the 
stability of the uuddsQ system within the chiral constituent quark model H used previously 
in the study of the H-particle [|7| and of heavy tetraquarks as well |1| . An essential difference 



with respect to Refs ||12| , pr5H l^] is that the spin-spin interaction of || is flavour-dependent. 
In Ref. [jn| we found that the Goldstone boson exchange interaction between quarks binds 



strongly both the ccqq system and the bbqq system. Within conventional models based on 
one-gluon exchange, the ccqq was found unstable and bbqq stable |T6| . 

In section 2, we establish the basis states required by the internal symmetries of the 
system under discussion and, based on simple group theory arguments, we indicate the 
most important basis states for a given isospin I and total angular momentum J. In section 
3, we briefly describe the Hamiltonian of the chiral constituent quark model used in the 
calculations. In section 4, we present our results and the last section is devoted to a summary. 

II. BASIS STATES 

In principle the GBE interaction contains all pairs ij of particles. But the exchange of a 
heavy pseudoscalar meson, between a light and a heavy quark Q, can in practice be neglected 
|i~8| . For example, it was explicitly shown in Ref. ]19[ that the dominant contribution to the 



masses of C=+l charmed baryons is due to meson exchange between light quarks and that 
the exchange of D and D s mesons is negligible, of the order of few MeV. Inasmuch as the 
quark-quark interaction ||,|| is inverse proportional to the masses of the interacting quarks, 
the exchange of B or B s meson can further be neglected. Neglecting the contribution of 
D, D s , B and B s mesons, the GBE interaction in the uuddsQ system, with Q=c or b, reduces 
to the interaction between light quarks q=u, d or s. 

We adopt this point of view, i.e. we neglect the contribution of heavy meson-exchange. 
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Then similarly to the study of the H-particle or of the NN system [2^] it is useful to have a 



qualitative insight about the uuddsQ system by first considering a schematic quark-quark in- 
teraction which simplifies the GBE interaction of Ref. ||, by removing its radial dependence. 
The schematic interaction reads: 

V^-C^Xf.Xfa^, (1) 

i<j 

where Af (F=l,2,...,8) are the quark-flavour Gell-Mann matrices (with an implied summa- 
tion over F) and a are the spin matrices. The minus sign of the interaction (|l|) is related to 
the sign of the short-range part of the GBE interaction, crucial for the hyperfine splitting in 
baryon spectroscopy. This feature of the short-range part of the GBE interaction is clearly 
discussed at length by Glozman and Riska [[| . A typical order of magnitude for the constant 
C x is about 30 MeV. 

In order to calculate the expectation value of ([!]) for the remaining light pentaquark 
system q 5 we have to give a classification of its states. In the colour space this system is 
described by the state [221]c compatible with the colourless state [222]c of q 5 Q. Here and 
below [f] stands for the corresponding Young diagram in the colour (C), spin (S) or flavour 
(F) space. 

We first assume that u,d and s are identical. If the quarks are all in the ground state the 
orbital part of the wave function is symmetric. Then the only q 5 flavour-spin state allowed 
by the Pauli principle is [32]^s. By using inner product rules J2T[| one can find the flavour 



[f]p and spin [f] s symmetries compatible with [32]^. These are listed in Table 1, together 
with the corresponding isospin I and spin S associated with these states. We also give the 
total angular momentum J = S + Sq of the q 5 Q system. The last column reproduces the 
expectation value of ([I]) in units of C x . This has been calculated using the formula given 
in the Appendix A of Ref. [p0[| , containing the Casimir operators of SU(6)ps, SU(3)p and 
SU(2)s- The multiplicity of a given IJ state is consistent with Table 3 (Y=2/3) of Ref. 
| |i2| |. From Table 1, one can see that the most favourable candidate for stability (the most 
negative eigenvalue of ([]])) should have 1=0 (flavour symmetry [22 1]^) and J=0 or 1. In the 



numerical calculations given below, different flavour symmetries will be mixed by the GBE 
interaction ([!]). 

In the diagonalization procedure of the GBE Hamiltonian given below, we truncate the 
basis for various sectors IJ, by retaining only the lowest states. The mixture with the others 
is expected to be small. Note also that the flavour spin interaction of type (|I|) does not mix 
[f]s 7^ [f'}s- In numerical calculations, we therefore restrict the basis states to the following 
ones : 

I=Q, J=Q or 1 

|1 >= |[221] F [32] S >; |2 >= |[32] F [32] S > (2) 

1=1, J=Q or 1 

|1 >= |[311] F [32] S >; |2 >= | [32] ^ [32] s > (3) 

1=2, J=Q or 1 

|1 >= |[41] F [32] S >; |2 >= |[5] F [32] 5 > (4) 



III. HAMILTONIAN 



The Hamiltonian to be diagonalized is 19 : 



# = £^+£jr - + £ W^) + £ (5) 

with the linear confining interaction : 

V conf {r lj ) = -h' i .\]Cr lJ (6) 
and the spin-spin component of the GBE interaction in its S77 F (3) form : 

v x(f*j) = { £ K(^)Af Af 
{f=i 

+ £ V^XfXf + V^r^X* + V^r^X") 9i ■ (7) 



F=4 



with A = y 2/3 1, where 1 is the 3x3 unit matrix. The interaction (|?|) contains 7 = tt, K, 77 
and rf exchanges and the form of V^r^) is given explicitly in Ref. || as the sum of two 
distinct contributions : a Yukawa type potential containing the mass of the exchanged 
meson and a short-range contribution, of opposite sign, the role of which is crucial in baryon 
spectroscopy For a given meson 7, the meson exchange potential is : 



a 2 1 I P -M r ij A 
Vjifa) = ~P~~o~a 04- p a 3 exp(-a 2 (r-r ) 2 )}, (7 = tt, K, V ,rf) (8) 

AtT 3 AmilTLj ' Tij \J TT 

For the Hamiltonian (||)-(@), we use the parameters of Ref. ||. These are : 

9irq _ 9rjq _ 9Kq _ g g« 9rfq _ ^ 
4/T 47T 47T ' AlT 

r = 0.43 fm, a = 2.91 far 1 , C = 0.474 fm _2 ,m M = 340 MeV. 

fi n = 139 MeV, fx v = 547 MeV, = 958 MeV, p K = 495 MeV. (9) 

They provide a very satisfactory description of low-lying nonstrange baryons, extended to 
strange baryons in in a fully dynamical three-body calculations as well. The latter 
reference gives m s = 0.440 GeV. For the masses of the heavy quarks we take m c = 1.35 GeV, 
mj, — 4.66 GeV in agreement with Ref. [18[] where these masses are adjusted to reproduce 



the average mass M = (M + 3M*)/4 of M = D and B mesons respectively. Note that within 
the spirit of the model of Glozman and Riska there is no meson exchange between a quark 
and a antiquark. 



IV. RESULTS 

First we discuss the spin S=l/2 baryons needed to calculate the threshold energy. Us- 
ing the Hamiltonian (^)-(^), we have performed variational estimates with a general wave 
function of the form ip ~ exp [—(ax 2 + by 2 )] where x — r± — 7%, y — (^1 + r*2 — 2r*3)/v / 3. 
To the nucleon mass, only tt, 77 and 77' exchange contribute. To the mass of A or E, there 
is a contribution from K-exchange as well. In the case of heavy baryons A c , A b , E c , S fe , S c 
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(10) 



and we neglect any meson exchange between a light and a heavy quark, in the spirit 
of the above discussion. The results are presented in Table 2, where for N, A and S we 
made the simplification a=b. The expectation values lie typically about 50 MeV above the 
experimental value. The results for S could be improved by taking a^b and those for A& 
made more realistic by tuning the mass of m b . But, for the present purpose, as will be seen 
below, these estimates are quite satisfactory. 

For the hexaquarks discussed here, it is useful to introduce the following system of Jacobi 
coordinates (where i=l,2,...,5 are associated with light quarks and i=6 with the heavy one): 

x = r i — f. 2 
V = r 3 - r 4 

^(fi +f 2 -f 3 - f 4 ) 
i = + r 2 + f 3 + f 4 - 4f 5 ) 

+ r 2 + f 3 + f 4 + f 5 - 5r 6 ) 
-Rcm = (wf*i + rnf 2 + mr 3 + mr 4 + m s f 5 + mgr 6 ) / (4m + m s + mg) 

Moreover, in the kinetic term only, we use the average mass : 

m = (4m + m s )/5 (11) 

for all light quarks. 

By assuming a ground state variational wave function of the form: 

* = (-) 3 (V exp [—{x 2 + y 2 + z 2 + t 2 ) - V] (12) 

TV TV Z I 

the expectation value E of the spin-independent part of the Hamiltonian becomes: 

6 h 2 c 2 A 5 . Ti 2 c 2 

E ° = — + ^ — H J~J72~ 

m 2a z m otr 

4„.„ n~ . /i,2 3, [l6) 



z 



w 



+-C[10d— + 5W — (- + -)] 
5 V 7ra V 07r a o 



For the matrix elements of the spin-dependent part, the fractional parentage technique f2~I 
has been used. Details are given in Appendix A. 
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In Table 3, we present results from the diagonalization of the Hamiltonian for 
the bases (@)-(f|). The column-heading M represents the lowest expectation value obtained 
in that sector. The next column gives the lowest theoretical threshold compatible with a 
given IJ. We also indicate the next to the lowest threshold whenever it is close to the lowest 
one. Note that using the experimental masses, Table 2, last column, the two thresholds 
may interchange their position. The last column is the lowest eigenvalue from which the 
threshold mass Mt obtained from Table 2 has been subtracted. The lowest eigenvalue is 
the equilibrium value obtained by minimizing with respect to the variational parameters 
a and b of ([12]) • 111 each case it turns out that b at equilibrium is approximately equal 
to the value of b given in Table 2, associated to the heavier of the threshold baryons. At 
equilibrium, we also find that the off-diagonal matrix elements of the GBE interaction are 
typically one order of magnitude smaller than the diagonal ones so that the lowest state 
does not change much through the coupling to the next state. A typical change is of a few 
MeV. This also proves that the truncation of the bases as in (0)-(0) is safe. The smallest 
M — Mt corresponds to IJ=00 or 01, as expected from the discussion following Table 1. In 
all cases M — Mt is positive and very large which means that none of the considered system 
is stable against strong decays. Actually, there is a substantial amount of repulsion, similar 
to the case of the H-particle 0. Thus within the GBE model used here, the heavy compact 
hexaquark uuddsb is highly unstable, contrary to the findings of Ref. [|ll[ or [|l!|. In the 
latter reference the system with 1=0, J=2 is bound by 13.8 MeV for the most favourable 
choice of the model parameters. 

The amount of repulsion found depends, of course, on the approximations used, and in 
particular on the treatment of the kinetic energy. A better treatment, where instead of the 
mass average (|TTD the kinetic energy is expressed in terms of the average of the inverse of 
the reduced masses , may slightly decrease the kinetic contribution but certainly will not 
change the above conclusion. The expectation is that the incorporation of the D or B meson 
exchange will not change the conclusion either. But in cases where uuddsQ is found to have 



a mass close to the lowest threshold, as for example Ref [|12l , a proper treatment of the 



kinetic energy is much more important in drawing a conclusion about stability under strong 
interactions. In a shell model description where one assumes that all the quarks are in an s 



state, as in p2[ , one can make an estimate of the kinetic energy of a system of six identical 
quarks from which one can subtract the kinetic energy of two separate clusters of three 
quarks each. This gives 3/4huj JF], i.e. a positive contribution, which may counterbalance 
the binding found in Ref. . 

One can also raise the question whether or not by increasing the number of heavy quarks 
the stability would increase. In | 18fl , we also investigated the system qqqqQQ in the Glozman 
et al. model ||, by using a similar procedure. There the most favourable configuration has 
1=0, J=l. For Q=c we obtained M-M T = 0.523 GeV and for Q=b, M - M T = 0.515 GeV. 
In both cases M T corresponds to the lowest threshold qQQ + qqq where m(ccu) = 3.514 
GeV and m(bbu) = 10.066 GeV. Therefore, these systems are also unbound in a compact 
configuration. 



V. SUMMARY 

In a chiral constituent quark model which successfully describes the light, strange and the 
presently known charmed and b-baryons, we have calculated the mass M of the hexaquarks 
H c (uuddsc) and H b (uuddsb) for various IJ sectors and compared it to the mass M T of the 
corresponding lowest threshold. We found that the smallest M — Mt value is associated to 
the IJ=00 or 01 sector. The quantity M — Mt is always positive and of the order of few 
hundreds MeV. This indicates that H c and Hf, cannot exist as compact systems. 

However, the existence of a weakly bound, molecular-type heavy hexaquark system, 
like the deuteron, cannot be excluded. The GBE interaction || generates a long-range 
attraction due to its Yukawa-potential tail and in principle, it can also produce a medium- 
range attraction from correlated two-pseudoscalar meson exchange. It is certainly interesting 
to pursue investigations in this direction, in a dynamical approach as the resonating group or 
the generator coordinate method, by incorporating six-quark states with orbital excitations, 
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as in the NN case 
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APPENDIX A: 



The calculation of the matrix elements of the interaction potential (0) between five 



light quarks is based on the fractional parentage technique described in Ref. ||21|| . Through 
this technique each five-body matrix element reduces to a linear combination of two-body 
matrix elements of the pair, 4 and 5, of quarks. This is possible due to the fact that each 
part (orbital, spin, flavour, colour) of the wave function is written as a sum of products 
of the first three and of the pair (45) wave functions. These functions have a definite 
permutation symmetry [f] and [f] respectively, where [f'] = [2] or [11]. The coefficients of 
these linear combinations have been obtained from the isoscalar factors of the Clebsch- 
Gordan coefficients of S$ as calculated in Ref. W2 . 



The calculation of the spin-spin matrix elements is trivial. Below we give the flavour wave 
functions which we derived in the Rutherford- Young- Yamanouchi representation, where the 
pair 45 is either in a symmetric or an antisymmetric state, as mentioned above. We denote 
by p and q the row of the 5th and 4th particle in a Young tableau and by pq and p~q a 
symmetric and an antisymmetric state respectively. Then the uudds states required in these 
calculations are: 

|[5]11> = yflij[ 3 ](uud)4>[ 2 ](ds) + Tjl^[ 3 ](udd)<f>[2](us) 

+ \fli , [3\{uds)(j)] i 2]{ud) + yf^ip^idds)^]^) (Al) 
+ \fk^} (uus)4> [2 ]{dd) 



|[41]11> = -y'|V'g 1] o (wrfs)0[2]M + ^^ 21] (dsd)(/)[2](uu) 

+\fli ) [2i]( usu ) ( l ) m( uu ) 
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(A2) 



I [41] 12 > = ^J^[3](uud)0[ 2 ](ds) + yjV[3]Mrf)0[2](Ms) 

^[3](mis)0 [2 ](oM) 



| [41] 12 >= -W^V[3](^)0[ii](^) - J^[ 3 ](wdd)0[ii](tis) 



|[32]22> = |^[ 3 ](MM(i)0[2](rfs) - |^[3](M(i(i)0[2](Ms) 

-y^[3]H^[2](«rf) + ^[3](gMs)0[ 2 ](ww) 
+ 3^^[3]( MMS ) < ^[2]( cM ) 



I [32] 12 > = §Vf 21] (udu)0 [2] (ds) + |^f 21] (<M0 [2] M 

"^V'pi] (uds)(f)[2] (ud) + ^ip^(usu)(/>[2](dd) 



I [32] 12 >= ^ -itf 21] (dvd)</>[ U ](us) - ^ -i)[ 21] (uds)<p {11] (ud) 

| [311] 1 1 > = ip[ Ul ](uds)(f)[2](ud) 

| [311] 13 > = ^f 21] (udu)0 [2] (<te) - ^^f 21] (d«d)0[2]W 
~\lla^2i\ ( uds )<t>V\( ud ) + \/S^f2i]( rfsrf )^[2]( MM ) 

|[311]13> = -\il)^(udu)(j)[ii](ds) + |V ; f2i](d M d)0[ii](-us) 

-vT^g°(^)^[ii](^) 



|[311]23> = y|V[3](^d)0[n](ds) - ^/|V[3](wdd)0[ii](«s) 

+ \/l' l l J m(uds)(f) [11 ]{ud) 
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|[221]23> = ^f 2 i]W0[2](^) + ^f 21] (^)0[2](^ 

\uds)(f) [2] {ud) - |Vf 2 i] 

-|Vf 2 l] (MSM)0[ 2 ](dd) 



1 ^f 2 'lj (^ S )^[2]( Mrf ) ~ \^[2l](dsd)(j)[2](uu) ( A12 ) 



|[221]23> = -f^y^l^mi^-f^^^^mttts) 



"IV'pi] (W^S) </»[n] (ltd) 



(A13) 



|[221]12> = ^ {lll] {uds)(j) { n]{ud) (A14) 

The states ^[2](«fr) and "?/>[ii](a&) are the symmetric and antisymmetric two-particle states. 
The ipi3](abc) is the symmetric three particle states. For mixed symmetry states ip^i] some 
care should be taken. As usually plj one has : 



if)pi] {udu) = -(udu — duu) (A15) 



1 
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ip^°(uds) = f — {2uds — 2dus + sdu — sud + usd — dsu) (A16) 



ip^f]°( u ds) = —-(usd + dsu — sdu — sud) (A17) 

However, it turns out that the states | [32] 12 > and | [32] 12 > contain the function ^pip 
the definition of which is : 

r [n (uds) = ^r { i{ (Ms) 



i. 



e. a linear combination of ( |A16| ) and ( |A17| ). While the states ( |A15| )- ( |A17| ) have a definite 



isospin the state ( |A18| ) is a mixture of 1=0 and 1=1. Thus the states | [32] 12 > and | [32] 12 > 



do not have a definite isospin. Therefore one has to project into a specific value of I in the 
calculation of matrix element of these two states. Calculation with or without projection 
indicate a difference of few MeV which is insignificant in the context of the present study. 
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TABLES 

TABLE I. Expectation value of the operator (||) in units of C x for all flavour [f]p and spin 



[f]s symmetries compatibles with [32] fs- The corresponding spin S and isospin I, together with 
the total angular momentum J are also given. 
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TABLE II. Variational solution of the Hamiltonian ©-(0) for spin S=l/2 low lying baryons 
compared to experimental masses 



Baryon 


Variatic 


jnal parameters (fm) 


Expectation value 


Experimental mass (GeV) 




a 


b 


(GeV) 


[23| 


N 


0.4376 




0.9696 


0.940 


A 


0.4486 




1.1654 


1.1156 


S 


0.4625 




1.2354 


1.193 


A c 


0.4683 


0.7099 


2.3268 


2.2849 




0.6320 


0.7201 


2.4889 


2.452 




0.5705 


0.6967 


2.5494 


2.470 




0.4678 


0.6509 


5.6147 


5.641 


s b 


0.6292 


0.6623 


5.7775 


? 




0.5683 


0.6339 


5.83629 


? 
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TABLE III. Masses M and mass differences M — Mt (Mr - the threshold mass) of heavy 
hexaquark systems uuddsQ (Q=c or b) of given isospin I and total angular momentum J 



System I J M Threshold M - M T 

(GeV) (GeV) 

N + E c 0.625 

uuddsc 0,1 4.144 

A + A c 0.652 

£ + A c 0.742 

1 0,1 4.304 

N + E c 0.785 

2 0,1 4.496 S + S c 0.772 



N + E b 0.619 

uuddsb 0,1 7.425 

A + A 6 0.645 



S + A 6 0.736 

0,1 7.586 

N + E b 0.780 

0,1 7.780 S + 0.767 
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